Criticality in the two-dimensional random-bond Ising model 
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The two-dimensional (2D) random-bond Ising model has a novel multicritical point on the fer- 
romagnetic to paramagnetic phase boundary. This random phase transition is one of the simplest 
examples of a 2D critical point occurring at both finite temperatures and disorder strength. We 
study the associated critical properties, by mapping the random 2D Ising model onto a network 
model. The model closely resembles network models of quantum Hall plateau transitions, but has 
different symmetries. Numerical transfer matrix calculations enable us to obtain estimates for the 
critical exponents at the random Ising phase transition. The values are consistent with recent 
estimates obtained from high-temperature series. 
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I. INTRODUCTION 

Two-dimensional (2D) models have played a special 
role in the theory of phase transitions.^ In 1944 Onsager's 
exact solution of the 2D Ising model gave critical expo- 
nents that were simple rational numbers, although dif- 
ferent than Landau theory. In the 1970 renormalization 
group (RG) calculations revealed exponents varying con- 
tinuously below an upper critical dimension, illustrating 
the breakdown of Landau theory. But it was unclear why 
the 2D Ising exponents and those for other exactly sol- 
uble 2D models were rational numbers. This fact was 
explained by the remarkable development of conformal 
field theory in the 1980s. Under the assumption of con- 
formal invariance at criticality, it .was possible to analyze 
a large class of 2D critical points.El Moreover, a first step 
was made towards a full classification of all allowed 2D 
phase transitions. 

Many physically important 2D phase transitions oc- 
cur in systems with quenched disorder. An example of 
particular experimental interest is the transition betweea 
plateaus in the integer quantum Hall effect (IQHE).El 
This transition has been successfully studied numerically, 
but so far has eluded analytic treatments via either RG 
calculations, exact methods or conformal field theory. 
Given the general power of conformal field theory in 2D, 
it has been surprisingly unhelpful in understanding such 
random phase transitions. 

In this paper we analyze a nontrivial 2D random phase 
transition which occurs in the simplest of all models: the 
2D Ising model with random bonds. Our approach is nu- 
merical, an.d closely parallels earlier work on the IQHE 
transitionfj We first map the 2D randoia Ising model into 
a variant of the Chalker-Coddingtona network model, 
which describes noninteracting chiral fermions. The ran- 
dom Ising transition corresponds to a fermion localization 
transition. We extract critical exponents numerically by 



standard transfer matrix methods. The values are in rea- 
sonable ag-seement with those recently obtained by Singh 
and AdlerO via high-temperature series. Unfortunately, 
this random Ising transition has also eluded any analytic 
treatment. 

To be more specific, consider the Ising model on a 2D 
square lattice, 



Hj — S ' JijSiSj, 



(1.1) 



with nearest-neighbor interactions Jy taken as random 
variables with a distribution Pj(Jij). For the simple dis- 
tribution Pj(Jij) = pS(Jij + J) + (1 —p)S(Jij — J), corre- 
sponding to a fraction p of antiferromagnetic bonds_the 
phase diagram, as established by various methods ,[Tlij is 
shown schematically in Fig. 1. 
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FIG. 1. Schematic phase diagram of the 2D ±J ran- 
dom-bond Ising model. 
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For small p there is a phase boundary separating 
the ferromagnetically ordered phase at low temperatures 
from the paramagnet. A larger p destroys the ferromag- 
netic phase, replacing it in high dimensions (d > 3) by 
a spin glass phase. In this case a multicritical point is 
expected at the coexistence point of all three phases. 
In 2D the spin glass phase is not present at T ^ 0, 
being destroyed by thermal fluctuations However, the 
multicriticaJ -point still exists on the ferro-to-para phase 
boundary.EJiij The multicritical point is unstable along 
the phase boundary, with RG flows as depicted in Fig. 
1, consistent with the (marginal) irrelevance of weak 
disorder at the purepSF) Ising critical point (p = 0). 
LeDoussal and HarrisE-3 have argued that the Mishimori 
line, along which the internal energy is analyticJLj passes 
through the multicritical point and coincides with one 
of the two RG scaling axes. The other scaling axis is 
tangent to the ferro-to-para phase boundary. Recently, 
Singh and Adlertl have obtained estimates for the two 
associated critical exponents, using a high-temperature 
series method. From a general point of view, this multi- 
critical point is of interest, being probably the simplest 
2D critical point which occurs at both finite temperature 
and finite disorder strength. 

Our paper is organized as follows. In Sees. II and III, 
we show that the random-bond 2D Ising model can be 
mapped, using a fermionic representation, to a variant of 
the Chalker-Coddington network modelo This mapping 
reveals a close similarity between the random Ising tran- 
sition and the IQHE plateau transition. However, due to 
a symmetry difference, the two transitions are not in the 
same universality class. In Sec. IV, we employ a trans- 
fer matrix approach to analyze the network model, and 
obtain estimates for the exponents at the random Ising 
multicritical point. Sec. V gives a brief summary and 
conclusion. 



II. FERMIONIC REPRESENTATION 

It is well known that the critical properties of the pure 
2D Ising transition are equivalent to a massless Majorana 
(real) fermion field JJ-Lp his studies of the bond-diluted 
Ising model, ShankartH constructed a model in terms of 
conventional (Dirac) fermions, by adjoining two identi- 
cal copies. For the random-bond Ising model, we show 
below that this procedure leads to a model of 2D chiral 
fermions, with a hopping_matrix element of random sign. 

Following ShankarjUl we consider a spatially 
anisotropic Ising model, retaining a lattice in one di- 
rection, but taking the continuum limit in the other 
(the "imaginary time" direction). The partition func- 
tion, when expressed in terms of a transfer matrix, can 
then be written Z =Tr exp(— /3Hid), where Hid is a ID 
quantum Hamiltonian and (3 is the system size in the 
"time" direction. The appropriate ID Hamiltonian for 
the pure Ising model is, 



Hid = Y}^< + K *«+il 

n 

where n label sites of a ID lattice and o~ a with a = x, y, z 
are Pauli matrices. This model exhibits a phase transi- 
tion when K\ — K2, which is in the (pure) Ising univer- 
sality class, as verified below. The transition also follows 
from a duality symmetry which exchanges high- and low- 
temperature phases. tH With the definition er^ = ^nf-n+i 
and <in = Y[j n <n fmi the Hamiltonian can be written in 
the form (2.1) with a a — > /i a and k\ <-*■ k 2 . 

The partition function is also invariant under Ki — > — Ki 
(i = 1,2), as seen by a spin rotation, cr" — > — er", with 
a = x,y for n odd and a — x, z for n even, which restores 
the Hamiltonian to its original form. This transforma- 
tion is equivalent to J — ► — J in in Eq. (1.1), mapping 
from a ferromagnetic to antiferromagnetic spin model. 
Thus randomness in the sign of the exchange interaction 
in Eq. (1.1) is equivalent to randomness in the sign of 
Ki. To describe the ±J Ising model, we thus consider 
interactions which are random functions of both n and 
imaginary time r, Ki — > Ki(n, r), taking either sign. 

A fermionic representation can be obtained by intro- 
ducing the Majorana fields 

m (n) = A=U «, m{n) = n «> ( 2 - 2 ) 

" m<n * m<n 

which anticommute, {^(n), rij(m)} = 5ijS nm . The 
Hamiltonian is quadratic when expressed in terms of 
these new variables, 

Hid = (-2i) }}Kir)i(n)r)2(n) - K 2 rii{ n ) r l2{ n + !)]• 

n 

(2.3) 

An identical copy of the system is constructed by in- 
troducing a new set of Majorana fields The Hamilto- 
nian obtained by summing the two, H = -^[HiDiVi) + 
HiD^i)], can be expressed in terms of standard (Dirac) 
fermion operators 

^i = ^={Vi+i^i), $ = _(7fc-i6) (2.4) 

as 

H = ^(-i«i)[^j(nh/>2(n) - ipl(n)ij}i(n)] 

n 

+(m 2 )[V4(n)V>2(n + 1) - ip\{n + l)^i(n)]. (2.5) 

Notice that H has a conserved U(l) charge: + w<2^2- 
At the pure Ising transition, there are gapless excita- 
tions in this conserved charge. With disorder present, 
the transition corresponds to a localization transition of 
these conserved fermions. 

To complete the mapping, we express the partition 
function as a path integral over Grassmann fields, 
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Z = J X>V^exp(-5), (2.6) 
where S is the Euclidean Action for H: 

Jt n 

(2.7) 

Reinterpreting imaginary time as a spatial coordinate, 
r — > x, the action S can be viewed as a 2D Hamilto- 
nian of chiral fermions, denoted Hid- To bring it into 
a canonical form, we define new right- and left-moving 
fermion fields 



energy eigenvalues of the 2D Hamiltonian. If a wave of 
energy E = is incident in the x direction, it will decay 
as exp(— |A|x), since p x = iA is pure imaginary. The 
decay length £ ~ corresponds to the correlation 

length of the pure 2D Ising model. The critical exponent 
is Vpure = 1 as expected. 

With disorder, the tunneling amplitudes, Kj, become 
random functions of position, n and x, and momentum is 
not a good quantum number. Nevertheless, at the Ising 
multicritical point (see Fig. 1), one expects the E = 
states of H2D to be extended, corresponding to an infinite 
correlation length. Away from criticality, one anticipates 
localized electronic states at E = 0, rather than a gap as 
in the pure case. 



TpRn = (-l)'Vl(n), TpLn = (-1)"^), 



(2.8) 



1>L = K- 1 ) n Mnl A 



-i(-l) n fa(n). (2.9) 
In terms of these the action becomes 



H 2 d = J dx ^2 bl>Rn ( i9 x)tpi 



(2.10) 



This Hamiltonian has a simple pictorial representation 
in terms of ID right- and left-moving fermion fields, cou- 
pled together by hopping strengths k\ and k^, as depicted 
in Fig. 2. The model closely resembles an anisotropic ver- 
sion of the Chalker-Coddington network model.Q In the 
next section we describe a lattice version, appropriate for 
numerical simulations. 
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FIG. 2. pictorial representation of the chiral fermion 
Hamiltonian (2.10). Neighboring chiral modes, propagat- 
ing in the direction of the arrows, are coupled via tunnel- 
ing matrix elements Ki. 

In the absence of disorder, with n\ and K2 constant, 
H2D can be easily diagonalized by transforming to mo- 
mentum space. The energy eigenvalues satisfy 



e 2 =pI 



2k\K2 cosp, 



(2.11) 



with p x the x component of momentum and p a trans- 
verse momentum in the range — 7r to n. The energy 
is a minimum when p — tt and p x — 0, and given by 
l-Elmin = =b I A I , with A = «i — K2- The pure Ising critical 
point occurs when A = 0. For A ^ 0, there are no zero 



III. NETWORK MODEL 

In pioneering work, Chalker and CoddingtonS intro- 
duced a network model to study numerically the transi- 
tion between IQHE plateaus. This model is essentially 
a lattice version of a chiral fermion Hamiltonian, similar 
to Eq. (2.10). The model consists of links and nodes, as 
depicted in Fig. 3. On each link there is a complex ampli- 
tude representing the fermion (electron) wave function. 
At the nodes, two incoming wave functions scatter into 
two outgoing ones, conserving probability. The nodes ana 
specified by an S matrix. In the original network model,cl 
the complex amplitudes acquired a random phase factor 
upon propagating along a given link, corresponding phys- 
ically to a (random) magnetic flux through plaquettes. 



jv„ v. 




FIG. 3. Representation of network model with width 
L = 4 and length iV = 3. The arrows indicate the direc- 
tion of wave propagation along links. The parameters 6i 
specify scattering at the nodes. 

On physical grounds, it is clear that a very similar net- 
work model should suffice for describing the propagation 
of E = waves of the Hamiltonian (2.10). In the appro- 
priate network model, the tunneling amplitudes are 
replaced by node parameters. Since we are interested in 
E = 0, there are no phase factors associated with the 
links themselves. The network model is specified by a 
transfer matrix T, taken, say, in the horizontal direction 
in Fig. 3. This matrix is decomposed into a product of 
matrices Mj , representing columns of the network, with j 
running from 1 to N, the length of the network. Each ma- 
trix Mj is a product of two matrices Mj = Aj(9i)Bj(02), 
representing two adjacent nodes in Fig. 3. The two node 
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parameters 9\ and Qi correspond to the two hopping co- 
efficients K\ and K2- 

The matrix representing tunneling at a given node is 
constructed to conserve the current or, equivalently. the 
U(l) charge. Following Chalker and Coddingtonjj the 
node in the dotted box in Fig. 3 is written as 

( w out \ _ ( cosh di sinh6>i \ ( v in \ 
\ Wi u J \ sm h 9i cosh 9i J \ v out J ' 

By construction, this matrix conserves the current, 
I Win | 2 + | Win | 2 = | w out | 2 + | Wout | 2 - Moreover, an incident 
wave, say, w- ln , is backscattered into u> ou t with probabil- 
ity tanh 2 (6*i ) . Since this tunneling probability is propor- 
tional to k\ in the 2D Hamiltonian (2.10), we make the 
identification 



tanh(6>i) <-> k^. (3-1) 

Although the Hamiltonian (2.10) is intrinsically 
anisotropic, the network model can be made 
invariantcl under a 7r/2 spatial rotation by choosing 
sinh(6>i) sinh(# 2 ) = 1, at successive nodes. 

It is instructive to briefly consider the pure network 
model, with constant node parameters. In this case, the 
transfer matrix can be diagonalized in momentum space. 
For a network with length N and width L, the total trans- 
fer matrix is T — Mq , with Mo an L x L matrix. Due 
to translational invariance, the eigenvalues of Mq, and 
hence T, can be labeled by a transverse momentum p. 
We denote these as A(p). For a given transverse mo- 
mentum p, these follow as eigenvalues of a simple 2 by 2 
matrix: 



( cosh 6*i cosh 02 + sinh 9\ sinh 02& %v smh 9\ cosh 02 + cosh 9\ sinh 02 e w 
\^ sinh 0\ cosh 02 + cosh#i sinh#2e -lp cosh£q cosh 02 + sinh#i sinh^e - ^ 



One finds, 



A±(p) =A±VA*-1, 



(3.2) 



with 



A = cosh(6»i - 2 ) + sinh^O sinh(6» 2 )[l + cos(p)]. (3.3) 

As before, at the pure Ising critical point a nondecay- 
ing mode is expected. Since the eigenvalues of the total 
transfer matrix are X N , this is only possible if |A| 2 = 1. 
This requires p — it and 0\ = 02, the expected condi- 
tion for Ising criticality. Specializing to p — n and the 
isotropic case, with sinh(6>i) sinh(#2) = 1, the eigenvalues 
take the simple form 



1 



(3.4) 



with A measuring the "distance" to the critical point: 



A = tanh^i) - tanh(0 2 )- 



(3.5) 



These eigenvalues describe the slowest decay of T ~ X N . 
The Ising correlation length follows as £ = l/ln(\ + ), 
and as expected varies as £ ~ 1/|A| upon approaching 
the critical point, A — > 0. 

We now incorporate randomness. Due to the identifi- 
cation (3.1), a change in sign of the Ising exchange, cor- 
responds to a sign change of a node parameter 0i in the 
network model. Thus the random-bond Ising model cor- 
responds to a network model in which the sign of the node 
parameters is random. To be specific, we choose the mag- 
nitude of the node parameters to be constants 0\ and 02 , 
satisfying the isotropy condition sinh(6*i) sinh^) = 1. 
The sign of at each node is chosen randomly, being 
negative with probability W and positive with probabil- 
ity 1 — W. The random network model we consider is thus 



characterized by two parameters: A, which measures the 
distance from criticality in the pure model, and W the 
disorder strength. 

This model differs in symmetry n from the original 
Chalker-Coddington network model,Q in which random 
fluxes were present through each plaquette, reflecting 
the breaking of time-reversal invariance by the magnetic 
field in the QHE. In the present case, there are no ran- 
dom phase factors. However, in the pure model with all 
node parameters positive, the fermion amplitude picks up 
a minus sign upon encircling any elementary plaquette, 
equivalent to a uniform flux 7r through each plaquette. 
Moreover, a sign change of a node parameter, changes 
the flux through two neighboring plaquettes by ir, as de- 
picted in Fig. 4. 




FIG. 4. A sign change in a node parameter 9 effec- 
tively changes by -k the flux penetrating the two neigh- 
boring plaquettes. 

Under the Ising duality transformation, A — > —A. 
Thus the A = line should correspond to the Ising fer- 
romagnetic to paramagnetic phase boundary, depicted in 
Fig. 1. Increasing the randomness W, with A = 0, corre- 
sponds to moving along this phase boundary away from 
the pure Ising critical point. Based on Fig. 1, we ex- 
pect the Ising multicritical point to correspond to some 
critical disorder strength W c . Numerical simulations, de- 
scribed in the next section, indeed support this scenario. 
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The scaling axis along the phase boundary corresponds 
to varying W with A = 0, whereas the Nishimori line 
corresponds to the line W = W c . 

For W > W c and A = 0, one expects the network 
model to remain critical, due to duality symmetry. How- 
ever, it is unclear what this regime corresponds to in 
the original lattice Ising model (1.1). It is conceivable 
that increasing W beyond W c (at A = 0) corresponds 
to moving along the low-temperature part of the phase 
boundary in Fig. 1, arriving at the T — fixed point at 
the maximum disorder strength, W max = 1/2. However, 
this interpretation is a bit problematic since W — 1/2 is 
"halfway" between the ferromagnet and antiferromagnet, 
and naively corresponds to p — 1/2 in Fig. 1. Perhaps 
the time continuum limit taken in Eq. (2.1), does not 
give a faithful representation of the lattice Ising model 
(1.1) for W > W c . 



IV. RESULTS 

The transfer matrix T of the random network model 
can be .-computed wurperically, following the work of 
ChalkerQ and others Jl3~E3 We have studied strips of width 
L ranging from L = 16 up to L = 128. Ensemble averag- 
ing is performed by taking very long strips, with length 
N up to 10 5 . Of the L eigenvalues of T/v, denoted Aj, L/2 
are greater than 1, corresponding to exponentially grow- 
ing solutions, and the others are less than 1, decaying to 
zero with increasing N. Due to the (statistical) parity of 
the system, these come in pairs, Xi = exp(±iV7j), where 
i = 1,2,..., L/2 and all the 7| are positive. Of interest 
is the smallest, 7 m i n , corresponding to the most slowly 
decaying mode. From this one extracts the correlation 
length as 

£l = — ■ (4.1) 

Tmin 

In order to extrapolate to the thermodynamic limit, 
L — > oo, it is convenient to consider the dimensionless 
ratio, 

Al(A, W) = (4.2) 

which is a function of the two control parameters W and 
A. Away from criticality (A ^ 0), £l=oo is finite, and 
so this ratio should grow and diverge as L — > oo. Repre- 
sentative data are shown in Fig. 7, where is plotted 
versus A for weak randomness W — 0.075, at various 
different system sizes. At the pure Ising critical point 
A = W — 0, the ratio A L is found to vanish identically 
even for finite L, due to the propagating mode with trans- 
verse momentum p = tt, described by Eq. (3.2). At a 
random critical point, one expects that Al will approach 
a finite constant in the thermodynamic limit, reflecting 
the infinite correlation length. 



In Fig. 5 we show data for versus disorder strength 
W, along the phase boundary, A = 0, at four different 
system widths L. Although there is significant variation 
with L, particularly for the smaller sizes, there appear 
to be two distinct regimes separated by a peak. For the 
largest width (L = 128) the peak occurs at a disorder 
strength W c w 0.08. For W < W c the ratio A L drops 
rapidly towards zero, the value at the pure Ising critical 
point. For strong disorder W > W c the ratio appears 
to be settling down towards a constant of order 1/2 at 
large L. This presumably corresponds to a strong disor- 
der critical point. Right at the critical disorder strength 
W = W c , Al is increasing slowly with L, and presum- 
ably eventually saturates. We thus identify the point 
(A = 0, W = W c ) with the random Ising multicritical 
point (see Fig. 1). 




FIG. 5. Dimensionless ratio Al vs disorder strength 
W along the phase boundary A = for four different 
system widths: L = 16(x), L = 32(A), L = 64(D), and 
L = 128(0)- 

This identification can be confirmed by extracting crit- 
ical exponents, ancLcomparing with the values obtained 
by Singh and AdlerB at the Ising multicritical point. Un- 
fortunately, precise values are difficult to extract due to 
the rather severe finite size effects, evident in Fig. 5. 
These are due (in part) to the small value of W c : With 
W c ~ 0.08, the typical nearest distance between two 
nodes with negative node parameters, 6i, is roughly 4 
times the network lattice spacing. 

Consider first the critical exponent corresponding to 
the scaling axis along the phase boundary, denoted v v . 
A natural finite-size scaling ansatz for A l takes the form, 

A L (A = Q,W) = /[L^^iW -Wc)}, (4.3) 

for large L and W — > W c . This form predicts that peaks 
in Al(W) around W c should sharpen up with a width 
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vanishing as SW ~ L~ 1 ^ Vp . The data in Fig. 5 are con- 
sistent with this trend, showing narrower peaks for larger 
L. To obtain a rough estimate for the exponent v p , we fit 
the peaks to a parabola. Denoting the curvatures of the 
parabolas as Rl, scaling predicts i?x ~ L 2 ^ p . In Fig. 6 
we plot log(i?x) versus log(L) and extract the exponent 
2/ Up as the slope of a fitted straight line. Fitting all four 
points gives an estimate v p ~ 2.2. However, the \ 2 P er 
data points decreases an order of magnitude if the small- 
est size L — 16 is excluded, which gives (dotted line in 
Fig. 6) v p ~ 2.45. Thus we estimate v p w 2.4, with a 
large error bar ±0.3. 
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FIG. 6. Curvatures Rl of parabolic fits to the peaks 
in Fig. 5, vs system width L on a doubly logarithmic plot. 
From scaling, the slope gives an estimate for the critical 
exponent 2/u p . 

The critical exponent along the Nishimori line can be 
extracted by sitting at W c and tuning A away from zero. 
In this case, finite-size scaling implies that 

Ax, (A, W = W c ) — F[LV V A], (4.4) 

The raw data for Ax, (A, W c ) versus A are shown in Fig. 
7. As expected, away from the multicritical point at 
A = 0, the ratio Ax grows with L, indicative of a finite 
correlation length. In Fig. 8 these data for A < 1/2 are 
replotted, rescaling the horizontal axis by L x l v . Based 
on the quality of the data collapse, we estimate i/«4/3 
with error bars ±0.1. 

Our estimates for the exponents compare favorably 
with those obtained by Singh and Adlera: v — 1.32±0.08 
and Vp k, (5/3)v « 2.2. This agreement gives one confi- 
dence that the peak in Fig. 5 does indeed correspond to 
the random Ising multicritical point. 

In addition to the multicritical point, we have tried 
to analyze the behavior at maximal disorder strength 
W = 1/2. As discussed earlier, it is unclear whether 
or not this point corresponds to the T = fixed point at 
p = p c in the random-bond Ising model (Fig. 1). Unfor- 
tunately, at W = 1/2 we are even more severely plagued 



by finite-size effects as A is varied. Specifically, the cor- 
relation length tends to remain very long, even well away 
from criticality, with A — > 1. In this limit, the fermions 
tend to become localized around plaquettes on one sub- 
lattice. But with W — 1/2, half of the plaquettes have 
zero flux, and can support (E = 0) states circling around 
them. For 1— A small, these states will be weakly coupled 
via tunneling, and may tend to percolate out to rather 
long scales. At this stage, we cannot conclude anything 
definite about the critical behavior of the network model 
at strong disorder. 
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FIG. 7. Dimensionless ratio Al = L/£l plotted vs A 
at disorder strength W = 0.075 for four different sys- 
tem sizes: L = 16(x), L = 32(A), L = 64(a), and 
L = 128(0) 

V. CONCLUSION 

To summarize, we have shown that the 2D random- 
bond Ising model can be fruitfully mapped onto a net- 
work model for chiral fermions. The network model is 
similar to the original Chalker-Coddington modelta used 
to study IQHE plateau transitions, but with different 
symmetry. Specifically, the model has node parameters 
with random signs, rather than random fluxes through 
plaquettes. The network model has been used to study 
the novel random multicritical point which exists on the 
Ising ferromagnetic to paramagnetic phase boundary. By 
implementing a numerical transfer matrix approach, es- 
timates for the associated critical exponents have been 
extracted, which pare consistent with those obtained by 
Singh and AdlerO from high-temperature series. The 
critical exponents are quite different from those at the 
IQHE plateau transition, indicating different universality 
classes, not surprising in view of the symmetry differences 
between the two models. 
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FIG. 8. Scaling collapse of the data from Fig. 7 with 
A < 1/2 for four different exponent values: v = 1.1, 1.2, 
|, and 1.5. 



What are the prospects for an analytic treatment of the 
random Ising multicritical point? Being in two dimen- 
sions, one might hope that powerful constraints from con- 
formal invariance would be helpful. Analytic approaches 
to the IQHE plateau transition have been impeded by the 
absence of critical behavior in ensemble-averaged single- 
particle Green's functions. The situation might be sim- 
pler at the Ising multicritical point, though, since critical 
properties are probably present in average single-fermion 
correlators. 
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